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Motivated by recent interest in spin-triplet superconductors, we investigate the vortex lattice 
structures for this class of unconventional superconductors. We discuss how the order parameter 
symmetry can give rise to U(l) xU(l) symmetry in the same sense as in spinor condensates, making 
half-quantum vortices (HQVs) topologically stable. We then calculate the vortex lattice structure of 
HQVs, with particular attention on the roles of the crystalline lattice, the Zeeman coupling and 
Meissner screening, all absent in spinor condensates. Finally, we consider how spin. orbit coupling 
leads to a breakdown of the U(l)xU(l) symmetry in free energy and whether the HQV lattice 
survives this symmetry breaking. As examples, we examine simpler spin-triplet models proposed in 
the context of Naj;Co02 ■ yH2 and Bechgaard salts, as well as the better known and more complex 
model for Sr2Ru04. 



I. INTRODUCTION 

A half quantum vortices (HQV) with vorticity h/Ae, 
which is half that of usual Abrikosov vortex with vor- 
ticity $0 = h/2e, presents an exciting example of frac- 
tionalized topological defects. Quantization of collective 
topological defects provides clear cut access to the na- 
ture of the ground state. For instance, the vorticity h/2e 
of Abrikosov vortex in a type II superconductor clearly 
shows that the circulation associated with the vortex is 
that of charge 2e Cooper pairs. Since the vorticity in a 
condensate is determined by the requirement of single- 
valuedness of the order parameter describing the con- 
densate, "fractionalization of vorticity" is possible with 
a multi-component order parameter when different com- 
ponents are allowed to wind separately. For instance, 
in a triplet superconductor, the additional Cooper pair 
spin degree of freedom can be free to rotate in plan^3^ 
giving rise to additional U(l) symmetry and an associ- 
ated spin winding number; cases where vortex fraction- 
alization is due to the U(l)xU(l) symmetry of different 
physical origin has also been studiecP. Hence observation 
of fractionalization of vortices can serve as an indicator of 
the structure of the order parameter in a given conden- 
sate. In addition, the recent proposals predicting non- 
Abelian fractional statistics for the composite of a HQV 
and the Majorana fermions bound at its core in the chiral 
triplet superconductors brought in recent rise in the at- 
tention and inte rest to the possibility of HQVs in triplet 
superconductor^^I^IMl This type of non-Abelian statis- 
tics was first studied for quasiholes in the spin-polarized 
f = 5/2 quantum Hall statc^ Therefore, if we want 
to obtain the same statistics for vortices in a spinful su- 
perconductor, the vortices should be HQVs so that there 
would be phase winding only for a single component. 

However, while there are a number of candidate 
triplet super condu ctors such as the single layer ruthen- 
ate Sr2Ru041ili2l, the cobaltate Na^^CoOs • yHjOi^ and 
organic superconductor^^, HQVs have never been ob- 
served in bulk systems, in line with the energetic stability 



issues raised by two of us in Ref. It was pointed out 
in Ref. [3] due to the absence of screening for the spin su- 
percurrent circulation required for HQV in triplet super- 
conductors, HQVs can be energetically unstable in bulk 
samples towards combining into full Abrikosov vortices 
despite their advantage in magnetic energy. Related con- 
siderations have appeared in the context of spin-triplet 
superconductivity in UPta by ZhitomirskjEl 

The main motivation of this work is to investigate the 
possibility of using high enough fields to generate a HQV 
lattice in triplet superconductors where the vortex lat- 
tice serves two purposes at once: 1) stabilizing HQVs 
at finite separation 2) providing an unambiguous signa- 
ture of its formation (halving of the vortex lattice unit 
cell). Experiments have already determined the v ortex 
lattice structure successfully at low fields in Sr2RuOj ^^ l ^^ l 
and the observed square lattice geometry was consistent 
with the theoretical prediction by one of the present au- 
thors based on a chiral triplet order parameter in Ref. [T5] 
However, Ref. [TS] considered the limit of strong spin- 
orbit coupling, which leads to vortex lattices of full quan- 
tum Abrikosov vortices. Recently, measurements of the 
Knight shift for t he field along the c-axi^i^l as well as 
ARPES dataP^Ml indicate that that spin-orbit coupHng 
is perhaps not so strong. Therefore, in this paper, we 
extend the studies of Ref. [TS] to allow for weak spin- 
orbit coupHng, leading to the possibility of HQV lattices 
for fields along the c-axis. Additionally, the organic su- 
perconductor (TMTSF)2C104 naturally has weak spin- 
orbit coupling and Knight shift measurements provide 
evidence for a spin-triplet state at high magnetic fieldsP, 
which is precisely the situation we consider here. We also 
provide an analysis of this case and the closely related 
case for cobaltate spin-triplet superconductor^^. 

In this paper, we study the energetics of different vor- 
tex lattice configurations. The key additional physical 
ingredient is the U(l) spin-rotational invariance of the 
Cooper pairs that arises in a magnetic field. This gener- 
ically leads to two different species of fractional vor- 
tices whose fractional fluxes sum to $0. When stable. 



these fractional vortices form interlacing lattices analo- 
gous to vortex. antivortex lattice configurations proposed 
by Refs. Gabay and KapitulnilJ^^l^ Zhangl^ in the con- 
text of a two-dimensional (2D) superfluid and the con- 
figuration in two-component Bose condensates proposed 
by Ref. Mueller and Hd^. 

The rest of the paper consists of the following. In sec- 
tion [IT] we give a pedagogical introduction to the symme- 
try properties of triplet OP. In particular, we will show 
how U(l)xU(l) symmetry can arise in the OP of such 
systems. In section |III| we discuss the form of Gibbs 
free energy that is allowed by various symmetries in the 
problem when spin-orbit coupling is not included. In 
section |IV| we provide the general theoretical framework 
for the vortex lattice phases. In section |V] we show that 
the lowest Landau level solution often provides an ad- 
equate description and we discuss this solution for the 
lattice of HQV's. In section |VI| we consider the effect 
of U(l)xU(l) symmetry breaking driven by spin-orbit 
coupling. In section VII we lay out predictions for how 
to detect the proposed HQV lattice structures and we 
conclude with a summary and outlook in section |VIII[ 



where the function /(k) depends on the specific odd 
angular momentum channel. The key simplifying fea- 
ture is the the orbital dependence is described by a one- 
dimensional representation encoded by /(k). There has 
been suggestions that the cobaltate Naa;Co02 • WH2O 
has a spin-triplet pairing through /-wave channeP^l, al- 
though dat a from the Knight shift experiments remain 
controversiaPSmi. in this case a common choice is /(k) = 
kx{kl — 3ky). However, the precise form of /(k) is not 
needed for our results. For the organic superconductor 
(TMTSF)2C104, there is a strong case that the system 
becomes a triplet superconductor under sufficient field 
H >20kOe. In this case, there are many proposals for 
/(k). However, again, the specific form is not needed for 
our results. 

For chiral order parameter symmetry expected of the 
Sr2Ru04, with d in the basal plane and no spin-orbit 
coupling, the order parameter has four complex degrees 
of freedom: 



A(k)= ^/(fc^,fc,) 



<T = ± 



A||,CT 




(2.3) 



II. THE TRIPLET ORDER PARAMETER 

The order parameter of a tri plet s uperconductor takes 
a matrix form in the spin spac^^^^^Il 



A(k): 



A„(k) 
AiT(k) 



An(k) 
Aii(k) 



dj^ I ^ dij 



+ id 



v 



(2.1) 

where the spin quantization axis is along the z direction. 
The triplet pairing requires A|| = A|| and a set of three 
complex functions of k, namely ((ia;(k), (iy(k), ^^(k)), 
were introduced to parameterize the gap matrix. When 
the three functions are collectively represented using a 
vector notation, the "unit vector" d(k) represents the 
symmetry direction (zero projection direction) with re- 
spect to the rotation of Cooper pair spin. In the pres- 
ence of the sufficiently high field along the c-axis, Zeeman 
splitting between electrons with opposite spins prohibits 



pairing, leading to A 



n 



0. In the d- vector nota- 



tion, this implies that the d-vector lies in-plane (perpen- 
dicular to the applied field). In the rest of the paper, we 
assume that the field is sufficiently large so that this is 
the case. In the context of strontium ruthenate, our re- 
sults apply for the field along the c-axis (this is also true 
for the cobaltates when we are discussing spin-orbit cou- 
pling in hexagonal systems). For organic and cobaltate 
superconductors, our results apply for the field along any 
two- fold or higher symmetry axis of symmetry. 

For non-chiral triplet order parameter symmetry, 
which is expected of the cobaltate Naa;Co02 •?/H2C(^and 
organic superconductorJSl^ the spin pairing gap matrix 
takes the form 



A(k) = /(k) 







(2.2) 



where the function /, like the non-chiral case discussed 
above, depends on the specific odd angular momentum 
channel and = + iaky . This is equivalent to 

d(k) = A+d_|_ exp(irw/3j.) + A_d_ exp(— m(/3g.), (2.4) 

where n is an integer (n — 1 for p-wave, n = 3 for /- 
wave) and ip^^ is the azimuthal angle associated with a 
unit vector k in the 2D plane (assuming quasi-2D setting 
with the angular momentum along the c-axis: 1 = z). 
Although A_ = for homogeneous chiral superconduc- 
tor, we will show that A_ 7^ often plays an important 
role in describing the vortex lattice structure of chiral 
superconductor . 



Eqs.(^2.2| and (2.4 1 clearly shows that under these cir- 



cumstances, the order parameter symmetry takes the 
U{1) X U{1) form, which can allow for HQV's with /i/4e 
vorticity associated with tt orbital phase winding and tt 
d-vector winding. 



III. THE GIBBS FREE ENERGY 

In order to identify stable vortex type and the lattice 
structure itself, we start with the Gibbs free energy in- 
cluding all the terms allowed by symmetry up to quartic 
order. As usual, the quartic terms determine the vortex 
lattice structure. Due to additional spin degrees of free- 
dom, the full expression for the Gibbs free energy involves 
a number of additional terms compared to singlet super- 
conductor case and it is instructive to consider different 
contributions separately: 



/ — fmag + fo ^ + /z ^ + fiJ + fsO + fhoin fin ■ i^-^) 



(4) 



p(4) 



2 



where /„ 



/^TT—hH /Att is the magnetic energy (the 



field h is the sum of the external field H and the screening 
field) , the superscript (2) indicates terms quadratic in OP 
and (4) quartic in OP 



Other than the conventional quadratic term /q 



(2). 



|A,; 



(3.2) 



the remaining quadratic terms in Eq.(3.1l are conse 



quences of additional spin degree of freedom for the 
triplet superconductors. The Zeeman coupling term 



f(2) 
J Z 



'k/i(|A| 



lA 12 



(3.3) 



plays an important role for the HQV vortex lattice by 
introducing a slight spin-polarization. This slight spin- 
polarization gives rise to two phase transitions as in the 
case of the A1/A2 phase of ^H^HEH. The inhomo geneous 

(2) 

part of the quadratic terms are of the form 



' c.c. 



(3.4) 



where Di — V i + (27ri/<i>o)Ai. For these terms, we re- 
quire rotational invariance up to the lattice symmetry 
with respect to orbital degrees of freedom only, which 
means that we require invariance with respect to rotat- 
ing DiS and the orbital component of A^'s. The complex 
structure of /j„ can result in the condensate wave func- 
tion of a different form than that of conventional SC. 
(2) 

fso ^^"^ quadratic spin-orbit coupling term assuming 
the spin-orbit coupling to be small and is discussed in 
Sec. |VI| In the presence of spin-orbit coupling the free 
energy have to be invariant under the combined discrete 
rotation of the orbital and spin degree of freedom spe- 
cific for the given lattice symmetry. For lattices with 
orthogonal or tetragonal symmetry, spin-orbit coupling 
may reduce the symmetry of the Gibbs free energy and 
tend to suppress HQV formation by introducing a length 
scale beyond which the HQV's cannot exist (this length 
scale diverges as the spin-orbit coupling vanishes). This 
implies that the vortex lattice spacing must be less than 
this length scale for the HQV lattice to appear. How- 
ever, we show that for spin-triplet hexagonal materials 
(specifically the two-dimensional Fg and FT representa- 
tions in the notation of Sigrist and UedaPS), even large 
spin-orbit coupling still allows for the existence of a frac- 
tional vortex lattice. This consideration may apply to 
Na,Co02-yH20. 

Among the quartic terms, represents the usual 

set of homogeneous terms. As is shown in Appendix \K\ 
certain quartic terms vanish in the weak-coupling theory. 
The terms that vanish are those that lift the energy de- 
generacy between the full quantum vortex (QV) and the 
HQV lattice. For this reason, we also include the inho- 
mogeneous quartic term . This term accounts for the 
difference between the spin phase stiffness p^p and the 
overall phase stiffness ps- Not only does this difference 



play an important role in the stability of isolated HQV's 
as it was shown in Ref . i3j it plays the role of tuning pa- 
rameter for the vortex lattice structure. 

With multiple systems in mind, we consider contribu- 
tions to the Gibbs free energy specific for non-chiral and 
chiral superconductors respectively. 



A. Non-chiral Triplet Superconductor 

As mentioned earlier, here we assume the orbital de- 
pendence of the gap function to be the same for all spin- 
triplet components. Formally, this means that the or- 
bital degree of freedom belongs to a one-dimensional ir- 
reducible representation of the point group. We apply 
our analysis to materials that have orthorhombic, tetrag- 
onal, or hexagonal point groups. One relevant example is 
a non-chiral triplet /-wave superconductor with hexag- 
onal symmetry, which has been proposed in the con- 
text of the the cobaltates Naa:Co02 • 2/H2O; in this case 



/(k) = K{ki - 3fc^) from Eq.(l2^. When the d-vector 
lies in the xy plane, the inclusion of spin-orbit coupling 
implies that formally this order parameter belongs to the 
Fg representation of the hexagonal point group (the con- 
sequences of spin-orbit coupling for this representation is 
discussed in more detail in Section VI I . 



With the in-plane spin rotational invariance, the rele- 
vant free energy within the assumptions stated above is 
given by 



fl^^ E i^.(lAAnP + |AAj,p), 



t—x,y,z 



/fl-/3i(El^^H' + /^2|A^nAj,P, 

i 

/(;:)= 7[AtA,,(D^A^).(D^Aj^)*+c.c], 



(3.5) 
(3.6) 
(3.7) 



For tetragonal and hexagonal point groups = Ky 
while for orthorhombic point groups, 7^ Ky. For the 
high field limit we are considering, it is possible to re- 
scale lengths in two directions perpendicular to applied 
field such that Ki = Kj for i ^ j (where Ki refers to 
the new coefficient in the re-scaled coordinates) . We will 
therefore ignore the difference between the Ki and as- 
sume that for orthorhombic point groups we are work- 
ing in re-scaled co-ordinates. The term f^^^ is not the 
most general such term allowed by symmetry. However, 
it is this term that allows the GL theory to give the 
same physics as in Ref. f31 Indeed, one can gain more 
insight into the vortex lattice solutions that minimizes 
Eq.(3.2l and Eqs.( 3.6|3.7 l by relating the coefficient of 
the inhomogeneous quartic term 7 to the stiffness ratio 
Psp < Ps which controls the energetic stability of a pair 
of HQV's^. Within the London approximation, the gra- 
dient terms in Eq.(3.5l and Eq.(3.7) amounts to phase 
bending energy which will be proportional to {ps + Psp) 
and (ps — pgp) respectively for A^ and |Ajj^|. Combin- 



ing Eq.(3.5l and Eq.(3.7l with the homogeneous solution 



3 



q;/(/?i — f32), we obtain the following 



relation between 7 and psp/ps'- 
Kiipi-h) 1 



7 



Psp/Ps 



a 1 + Psp/Ps 



(3.8) 



Hence 7 > would imply stability of HQV's and double 
transitions into two possible vortex phases: a lattice of 
ordinary Abrikosov vortices and a lattice of HQV's. This 



transition is determined by the (i^ term of Eq. (3.6 1 and 
the 7 term of Eq. (l3j|. 



B. Chiral Triplet Superconductor 

With the ruthenate Sr2Ru04 in mind, we consider a 
chiral triplet p-wave superconductor with square symme- 
try for which /(fco-) — 



A(k) = ^{k^ + iaky) 



iaky in the Eq.(2.3l with 
A 











(3.9) 



where A^^^ (s =X{ , II and cr = ±) form expansion param- 
eters for the Gibbs free energy. In terms of the d-vector 
notation d = Mvxxkx + Vxvk-u 



■ yiVyxkx 



Vyyky), 









^Vyx 


-??ay)/2, 




{j}xx 


^Vxy 


^ ^Vyx 
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{j}xx 


^ ^^xy 


+ iriyx 


+ Vyy)/'2, 




{j}xx 


^Vxy 


+ iriyx 


-^yy)l2. (3.10) 



Formally, without spin-orbit coupling, this order pa- 
rameter is a direct product of a orbital represen- 
tation of the tetragonal point group and the in-plane 
vector representation for spin rotations. When spin- 
orbit is included the order parameter contains the four 
different one-dimensional representations of the tetrag- 
onal point group. In the case without spin-orbit cou- 
pling, the relevant free energy for this representation can 
be constructed using the known free energy for the 
representatioiP^, we list below Z-^^-*, and . Be- 

fore listing /j^/, we note that this free energy term should 
respect the C4 symmetry on the xy plane only for the or- 
bital degrees of freedom: 



(Z?„I?j,,A,.+,A,._) ^ (i?y,-D„iA,.+,-iA,._). 

(3.11) 

However, for simplicity, we consider cylindrical symme- 
try in the orbital degrees of freedom, this does not signif- 
icantly alter the arguments below. This symmetry gives 



f(2) 



J2s fin'^ where 



= ifi(|DA3,+ 



|DA,^-n 

{(i?,A,,_)(I?,A,,+)*-(i?,A,,_)(A,A,,+)*}/2 
i{{D,A,^^XDyA,,+y + {DyA,^^){D,A,,+y}/2 
i{{D,A,,+XDyAs,-r + iDyA,^+lD,A,^^)*}/2] 
K4\D,A,^+\^ + \D,A,^^\^). (3.12) 



In addition, the following term is also allowed by sym- 
metry 



sK^hy2{-\A 



s-l-l 



|A.-|^), 



(3.13) 



and it stabilizes this in-plane chiral phase for strong 
enough magnetic field (note the similarity to the Zee- 
man term for the condensate spin degrees of freedom). 
As for the homogeneous quartic terms. 



J hom 



^[/3i(|A,,+ |4 + |A,_n/2 + /3i|A,+ nA,_p] 



<T = ± 



/33[(Att,+Aj^, 



.)(AtT,-Aj^,+)*+c.c.]. (3.14) 



(32, and P3 terms originate from interaction between 
spin up-up pairs and down-down pairs. Again, for sim- 
plicity, we written the free energy in the limit of a cylin- 
drical Fermi surface. Lastly, we have 



(4) 



7 J2 [An,.A^^,.(Z)A„,.) . {DAl^„) + c.c 



+ i 



I][^TT.-Aii,-.(-DA„,„j 

<T = ± 



)-(r>AJ^,_J-Hc.c.]. 

(3.15) 



Note that the form of Eq. (3.151 is consistent with the 



form of the interaction terms in Eq. (3.141. Again, this 
is not the most general term allowed by symmetry, but 
it is the minimal term that captures the physics in the 
London limit described in Ref . El 



IV. DETERMINING THE VORTEX LATTICE 
STRUCTURE - GENERAL CONSIDERATIONS 

We consider the vortex lattice phases near the upper 
critical field to map out the stability condition for HQV 
lattice phases. As usual, the first step towards determin- 
ing the vortex lattice structure is to identify the eigen- 
states of the linearized Ginzburg-Landau (GL) equations. 
In order to obtain a linearized GL equation we take a 
variation of the quadratic terms in the free energy, for 
example: 



f(2) , .(2) , .(2) 
h) ^ iin ' J Z 



-a^lA^f -^M|AttP-|Aj^P) 
i 

[Kjk;im{DjAi){DkA„,y + c.c] 



(4.1) 

with respect to one component of the order parameter 
A*. This gives an equation of the form 



aA,=Kli.^^,DkDiK-kHA 



(4.2) 



(note that we are ignoring the difference between h and 
H in this approximation). Since the gradient terms in 



4 



Eq.(4.2l cannot in general be reduced into a + Dy 



form, the lowest Landau level wave functions are not suf- 
ficient for calculating the condensate wave function in 
general. However, the solution of this equation can still 
be expressed in terms of Landau level wave functions: 

</)„(r) = [2'Vi/2(„!)]-i/2 

X J2 'Zme''^-'"^'e-(^'-'=")'/2i/„(y' - fc„), 



V. LOWEST LANDAU LEVEL SOLUTION 

In the bulk of this section we provide detailed analy- 
sis of the lowest Landau level solution for the non-chiral 
triplet superconductors and briefly comment on the chi- 
ral case in subsection IV CI In this case the relevant free 
energy (for orthorhombic, tetragonal, and hexagonal ma- 
terials) is 



(4.3) 

where iJ„ is the Hermite polynomial of n*'' order, x' 
and y' are x, y coordinates in the unit of the mag- 
netic length / = (<I>o/27rif)^/^. This is because A's 
can be expressed as a linear combination of the raising 
and lowering operators of the Landau levels II-i-, since 
n± = l{D,±iDy)/V2. 

This wave function describes a vortex lattice when 
|(/)„(r)| is periodic in the lattice vectors ai ~ a/(l,0) and 
3.2 — bl{cos9,sm9), and (/)„(r) vanish at toiHi -|- 7712^1 
when mi and 7712 are integers. This requires 

km = 27r(m- l/2)/a= (777- l/2)\/27r(T 

iirm(?+l — m<;) 



(4.4) 

where a — (b/a) sinO and c; = (b/a) cos 6. Note that we 
used the flux quantization condition ab sin 9 — 27r. For 
a lattice of HQV's we need to consider a second lattice 
that is translated by It — l{Tx,Ty) with respect to the 
first lattice. For the wave function of this lattice, we can 
useP 



0n(r) 



.(r-x), 



(4.5) 



the phase factor being chosen so that n_0o(r) = 0- This 
latter condition ensures that the translated eigenstates 
have the same gauge as the untranslated eigenstates. 

The formalism considered here gives us not only the 
energy due to interaction between vortices but also the 
core energy of vortices as well. This is because our vortex 
lattice wave function gives full description of the core 
regions. From the linearized GL equation we use here, 
a full quantum vortex is merely two HQV's of opposite 
spins coinciding at a same point. This means that the full 
quantum vortex core energy, if we ignore the cross term 

between two spin components in fll^J^.^ approximately 
twice the core energy of a HQV. If the HQV core energy 
is actually larger than this, that would make stabilization 
of the HQV more difficult, i.e., the largest allowed value 
for Psp/ps for the HQV lattice would be smaller than 
what we obtain through the formalism used here. 

The lattice structure can be determined by flnding 
((T, c;, t) that minimize the free energy expectation value. 
For this, we first set the amplitude of OP to minimize 
the energy for given (cr, <j, t) (the amplitude depends on 
Hc2 — H), and then compare energy for different values 
of (cr, r). To determine these structures we will need 
to take the spatial integral of the product of four Landau 
level wavefunctions. We have computed these integrals 
in the Appendix [Cl 



/= E 

s=TT,i4. 



-alA, 



7[A* Ai^(DATT)-(DATT)*-fc.c] 



t—x.y,z 
2\ 



Hh 

Stt 47r ' 



(5.1) 



As mentioned before, we assume that we have re-scaled 
lengths so that we can take = = Ky = K. First 
look at the upper critical field problem. The linearized 
GL equation is: 

;2 / TTi2 \ 

^A,= h+2H+H_-,s — «j A,. (5.2) 

The largest Hc2 occurs when the A^ are in the lowest 
Landau level. This leads to two possible values for the 
upper critical field, 



H. 



a<i>o 



c2 



2T:{K±HPk) 



(5.3) 



We assume that k > so that A|| has the larger Hc2- 
We do not assume that the splitting between these two 
critical fields is large since it is given by the small Zeeman 
term k. 

Now consider the expectation value of / in terms of the 
Landau level wave functions. Applying Eq. (4.2 1 to the 



first two terms of Eq. (5.1 1 gives the lowest Landau level 
solutions A|| = C^cjiQ and A^ = Ci 



/jj^e*^"(^o- Inserting 
this solution as determined dX H — Hc2 gives 



-^(j.5A)-H/3i(|(/.oh(C^ + q\)/2 



{hi 



IP 

8^ 



1 - 



K + H^2l k 
K - Hc2pR 



p2 



(5.4) 



where j is the supercurrent, bh. is deviation of the vector 
potential from what we would have for the h = i?c2, and 
/is is the screening field of superconductor. Note that in 
this approximation j is calculated solely from quadratic 
terms, ignoring 7 terms, and by the Maxwell equation 
W X hg — iirj /c. More specifically 



—ckV X z(\A 



11) + c.c] 



TTl 



lAjiP) 



(5.5) 



5 



Since V x SA = z{hs + H — Hc2) and, from Maxwell's 
equations, W x hg — inj /c, partial integration leads to'^ 



1 



1 

47r 

(M) 

47r 



■{hs 

Hc2 



47r 



(5.6) 



Meanwhile, 
7, we set h = 

Hc2 



if) 



when we calculate the expectation value of 
= Hc2- This leads to the free energy of 

47r ^ " Stt Stt 



1 - 

i/2 



K 



K - Hr7pk 



p2 



f3 



TT 



2|^o|')-^</'o|'l0i|^)-/32^'/'o|'|^o|^ 



(/)• 



(5.7) 



The screening field hs can be calculated by assuming 
|A| cx {Hc2 — Hy/^ near the second order phase transi- 
tion aX H — Hc2- Here we will deal only with 0(|A|'*) 
(or equivalently, 0(1 — H/Hc2)^) and ignore higher order 
terms. This allows us to calculate j, and consequently 
/is, solely from quadratic terms. In the lowest Landau 
level, this yields the screening field: 



$0 



■—Auk 



Ctrl' 



47rK(| A|| 



111 

471^)0^1 00 1 



(5.8) 



Inserting Eq. (5.8 1 into Eq. (5.7 1, the free energy takes 



the following form: 

(/) = -diCl-d2Cl + l3iC^^ + P2Cl^+lizClCl (5.9) 

with terms quadratic or quartic in C|t and with co- 
efficients that are independent in general (we will further 
specify these coefficients is the next two subsections). In 
the absence of screening fields, Zeeman-fields and the 
term proportional to 7, the form of the free energy in 
Eq. (5. 7 1 is similar to that examined in Ref. 2A. in the 



context of two-component Bose condensates (spin-half 
spinor condensate) .In that case, the vortex lattice struc- 
ture is solely determined by the competition between the 
Pi and (32 terms of Eq. (5.1 1; the /3i term determines the 



interaction energy within each vortex lattices, and the (ii 
term the interaction energy between two fractional vor- 
tex species each forming lattices. Specifically, the quar- 
tic term — /32(|0oPl0P)C'||C|| determines the stability of 
the HQV lattice, li (32 < 0, then a HQV lattice is the 
ground state. If (32 > 0, then full quantum vortex lat- 
tice is the ground state. In Appendix [X] we show that 
/32 = in weak-coupling theories, so that the two lattice 
structures are degenerate. In the rest of the paper, we 
will focus on aspects that are unique to triplet supercon- 
ductors in subsections IV Al and IV Bl the effects of the 
screening fields fl^ 1 the Zeeman- field. 
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FIG. 1: Real space vortex lattice structure 



A. The effects of screening and 

Here we look into the effect of screening. Ignoring the 
Zeeman field (in weak-coupling theories, the Zeeman field 
is vanishing in the clean limit), we have symmetry be- 
tween C|| and Gil, giving us 



and 



ai — a2 



(3i = f32 = (3 



-{m\ ) 



$0 



P _ 8^3X2 

2 



(5.10) 



(5.11) 



which means that the free energy in Eq.(5.9 1 takes a sim- 
pler form: 

(/) = -d{Cl + Cl)+ p{C^^ + Cl)+ P^C^Cl, (5.12) 
with 

/?3 = ^(i<^op(i0op-i0in)-(^^+/?2) (i^on^op) 

Now the free energy of Eq.(5.7l can be minimized by 



choosing = C^^ = d/{2l3 + P3) (which gives |A| cx 

{Hc2 ~ 11)^1^ as mentioned), giving us the free energy 
expectation value 



(/) = 



8^ 2^ + /33 



(5.14) 



Eqs. (5.131 and (5.141 allows for understanding the role 



of both screening and j'^^ . The screening affect the vor- 
tex lattice structure throug h the dependence of terms 
proportional to in Eqs. (5.111 and (5.13) on the lat- 



tice structure parameters <j, a and r. Since all quartic 
expectation values depend on the lattice structure, the 
lattice structure will be determined through minimizing 
(2/3 -I- /Ja). Since the magnitude of term in Eq. (5.131 
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ItJ/laj + ajl 



1/2 
1/3 i 








Ia2l/lail 



0.182 0.370 0.926 
-(|3^+167£k2/^^)/P, 

FIG. 2: Transition between different vortex lattice structures 
when 7 = 0. (a)-(e) label the vortex lattice structure shown 
in Fig. ([l}. 



is larger for full quantum vortices (Eq.(5.11) is not af- 
fected), screening tend to disfavor HQV lattices, in line 
with earlier observation for isolated HQV'^. However, 
since weak-coupling theories lie near the point (3^ — Q we 
expect the screening effect will put the physical system at 
a fine balance between interlacing lattices of HQV's and 
the ordinary Abrikosov vortex lattice. Hence it should 
be possible to observe the transition between the two 
phases upon small change of field and temperature. In- 

("4) 

deed, the interaction flJ plays this role. In particu- 
lar, we numerically find that for 7 > 0, then this term 
tends to favor HQV lattices. A positive sign of 7 occurs 
when Psp,Ps and this is to be expected in spin-triplet 
superconductors'^. Note, that unlike the other contri- 
butions in /33, the contribution from fj^^^ vanishes as 
T Tc- Consequently, this term can drive a field and 
temperature dependent transition between a HQV lat- 
tice and a full quantum lattice. Figures 1 and 2 show the 
phase diagram for 7 = 0. Note the similarity between the 
calculated phase diagram and that found in the context 
of two-component Bose condensateJ^SMI^ 



B. The effects of the Zeeman term 

For simplicity, we consider here the effect of Zeeman 
field alone ignoring screening and setting 7 = (ignoring 

fin)- presence of Zeeman field, the free energy 



terms: 



(/) = -diC2 - d2Cf^ + ^(C^ +Cf^) + /SsCf^Cfi (5.15) 



where ai — a+^ + HR, 0:2 = a+^ — Hk, (3 = (3i{\(f>o\'^) , 

and — — /32(|<?!'o|^|<^oPP)- The Zeeman field has two 
main consequences: (i) it typically leads to two phase 
transitions. In the first phase 7^ and = 
and in the second phase, both components are non-zero. 
The first phase is analogous to the ^He Ai phase, with 
a non-unitary spin-triplet order parameter. However, 
weak-coupling theories prefer unitary spin-triplet states 
and this drives the second transition, (ii) In the frac- 
tional vortex lattice phase where both components are 
non-zero, the magnetic flux contained by isolated frac- 
tional vortices is no longer a half-integral flux quanta. 
Instead, two types of vortices each carry fractional flux 
values of 



C2 



(5.16) 



The double transition is possible ii 2(3 + > (note 
again that weak-coupling theories yield — and /3 > 
0), there can be two transitions with a second transition 
appearing at a temperature 



c2 



Af3 K + Hc2fk 



2/3 + PsKi- Hc2pfi 



(5.17) 



Eq. (5.9 1 would have different coefficients for two quartic 



In the high temperature phase, the vortex lattice is 
hexagonal and, at the second transition, the lattice will 
remain hexagonal and the second component will either 
coincide with first or be displaced half a hexagonal vortex 
lattice vector from the first. As temperature is further 
reduced below the second transition, the lattice will con- 
tinuously deform, asymptotically approaching the phases 
presented in the subsection V A| (those shown in Fig. 2) . 
The resulting phase diagram is qualitatively shown in 
Fig. 3. 

Both consequences of the Zeeman field stem from 
breaking the additional Z2 symmetry that is present 
when fii — a2- In general, the existence of fractional 
vortices is the result of the C/(l) x U{1) symmetry of the 
free energy. When there is additional Z2 symmetry due 
to ai — 012, the fiux contained in each fractional vortices 
are restricted to be half the fiux quantum since the two 
components of the order parameter are no longer degen- 
erate in a magnetic field. In Section VII, we will see that 
this helps us distinguish a lattice of HQVs from a lattice 
of full quantum vortex. 



Chiral Triplet superconductors: lowest Landau 
level solution 



The chiral triplet superconductor with tetragonal sym- 
metry, because of the inhomogeneous quadratic terms we 
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H 




Normal vortex phase 

Fractional flux vortex phase 




FIG. 3: The role of the Zeeman field is to cause two transitions 
and to change the HQV lattice pheise to a a phase with two 
types of fractional vortices in which the two fractional fluxes 
sum to <l>o . Far below the second transition it is expected that 
these fractions will be well approximated by 1/2. 



have already seen in Eq. (|3.12l 



+ K2[{{D,^s,+lD,^s,-T~{Dy^s,+lDy^s,-)*}l'2 
+ {(Z?,A,,_)(i^,A,,+)*-(i?j,A,,_)(D^A,,+)*}/2 
+ j{(D,A,_)(Z?j,A,,+)* + pyA,,_)(i?,A,,+)*}/2 
- i{(D,A,,+)pyA,^_)* + (i?yA,,+)(Z?,A,,_)*}/2] 
+ Ki{\D,As,+ ? + \Dz^sA''), (5.18) 
has a much more comphcated quadratic free energy, 



Jo 



(2) _ 



E 

s=TT,i4- 



ha(|A,,+ p + |A,,_p) + /^^^)], (5.19) 



even when we exclude the Zeeman field and any spin- 
orbit coupling. 

Due to these inhomogenous quadratic terms, we cannot 
put both chirality components in the lowest Landau level. 
This is due to the presence of [{DxAs,a){DyAs^-a)* +c.c] 



terms in 



(2. 



The above quadratic free energy of 



Eq. (5.191, together with Eq. (3.131 that gives us the 



energy splitting between two chiralities, leads to the lin- 
earized GL equation 



A,_ 



is:i(i+2n+n_ 
K2I1I 



-5K K2III 



(5.20) 



Eq. 



However, a lowest Landau level solution that satisfies 
(5.20 1 may still have the highest Hc2 and therefore 
be possible^^. This would lead to a nonzero order pa- 
rameter for only one chirality - (As+,As_) = C{0,<j>o) - 

and requires SK < —jj^- In this case, the vortex ener- 
getics of the chiral triplet superconductor is identical to 



that of the nonchiral triplet superconductor, for inserting 
this lowest Landau level solution into the full Gibbs free 
energy leads us b ack t o Eq. (5.4). 

However, Eq. (5.20) can also give us a solution with 

and A, 



Landau level mixing; in this case, both As-|_ and As_ are 
nonzero. We will present discussion on this Landau level 
mixing in Appendix IS] 



VI. ROLE OF SPIN-ORBIT COUPLING 

A. Generic Case 

If the material has orthogonal or tetragonal symme- 
try (though not necessarily true for hexagonal symmetry, 
which is discussed in the next subsection) , then there will 
exist spin-orbit coupling terms of the type 



eA 



TT ^ii- 



(6.1) 



Such terms break the U{1) x U{1) symmetry and con- 
sequently isolated fractional flux vortices are no longer 
stable. Nevertheless, a fractional flux quantum vortex 
lattice can still exist, provided that the separation be- 
tween vortices is less than ^^o defined through = K/e. 

For completeness, we write here the spin-orbit coupling 
terms that appear in the context of a chiral spin-triplet 
superconductor. While we do not include these terms in 
calculations, they may be useful in other contexts. Due 
to the tetragonal C4 symmetry, homogeneous spin-orbit 
coupling terms should be invariant under the transfor- 
mation 

(A||,+, A|^ Ajj_^+, Aj_j_ _) 
(-^TT.+ '^TT -'^ll.+ '-^-U -)■ 
To the quadratic order, this condition is satisfied bjEll 



/^^ = ei(|AtT,+ P + |A^,-P 
+ e2[(A^,_)(Aj^,+)*+c.c. 
+ e3[(ATT,+)(Au,-)*+c.c. 



Att,-P-|Aj,,+ |2) 



(6.2) 



We note her e that e^'s can be estimated from the recent 
ARPES datcPSHa. 



B. Hexagonal Materials 

For hexagonal materials, there exist spin-triplet pair- 
ing states for which no such terms such as that in the 
above equation appear. These states belong to the two- 
dimensional representations labelled and Fg in the 
review article by Sigrist and UedaP^. Consequently, these 
materials need to be considered more carefully. 

We will now show that in hexagonal materials, a little 
away from Hc2, spin-orbit coupling does not break J7(l) x 
[/(I) symmetry. For hexagonal materials, the only term 
that exists in the GL free energy that is due to spin-orbit 
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coupling is (note that the inclusion of this term gives rise 
the complete free energy found that is found in Sigrist 
and Ueda for the F^g representations): 

fso = K,o{iD^.Au)iD,A^^r-{DyAu){DyA^^y 

+ (i^,An)(i?.Au)* - (i?,Au)(i?,An)* 
- z[{D,A^^)iDyA^^r + iDyA^^){D,A^^r] 

+ z[(i?,An)(i^,Au)* + (I?,An)(i?.Au)l}/2 

(6.3) 

With the field along the c-axis, the solution to the 
quadratic problem satisfies 



1< \ A 



Hi 



1 + 2N + K, 



where K, 



K 



and Kso — Kgo/ K. All the eigenstates 
for this problem can be found analyticallj^^Hl. Typically, 
\Kso\ « 1, so we will be interested in the eigenstates 
that contain the lowest Landau level (which will minimize 
the free energy when Kso = 0). The two relevant eigen- 
states that we wish to keep are: (A||, A^) = ((pQ, e(f)2) 
and (A||,A|.|) = (0, 0o) where e is proportional to 
Kso- Note that unlike in the subsection V C we can 
keep both solutions because while we examined H ~ Hc2 
in that subsection, we are a little away from Hc2 in 
this subsection. We therefore write A|| — 7i(/)o and 
^ii — li^4'2 + 7200 to include these two eigenstates. 
For simplicity, we ignore screening and the Zeeman field 
to find the following free energy 

(/> = -(1 - H/Ho2,l)\ll? - (1 - H/Ho2,2)\l2? 

+/3i[|7i|'(l'/'o|'> + (l7ie02+72<^o|'>] 
-/32(|7i0on7ie02 + 72^oP) (6.4) 

where ffc2,i (* = I72) is the upper critical field for eigen- 
state i. Since spin-orbit coupling is expected to be small, 
this implies that e << 1, so keeping to linear order in e 
yields: 

(/) = -(1 - i//i?c2,l)l7l|' - (1 - H/Ho2^2)\l2? 

+/3i(l0o|')(|7i|'+l72n+/32(|</'on0oP)l7in72|' 



-e[/3i|72|'727i*(l^o|'0o0;) + c.c] 
-e[/32|7i|'7i72(l'/'ol'020S>+c.c] 



(6.5) 



Without the last two terms, this theory is the same 
as that found for non-chiral spin-triplet superconductors 
with a Zeeman field but without any spin-orbit coupling. 
At the upper critical field, one of the two components 71 
or 72 order and the vortex lattice will be hexagonal (this 
conclusion is correct even when including terms that are 
second order in e). As the temperature or magnetic field 
is reduced, the last two terms in Eq. (6.51 can play an 



important role. These two terms break the U{1) x C/(l) 
symmetry of the theory and therefore will tend to remove 



any HQV lattice phases. However, the spatial averages 
{\4>o\'^ 4>q4>2) and (|(/'oP02</'o) vanish for a hexagonal vor- 
tex lattice (loosely speaking, this follows from noting that 
(pn picks up a factor e™*^ under a rotation about z and 
that a hexagonal vortex lattice is symmetric under rota- 
tions of 7r/3) . The hexagonal symmetry of the materials 
conspires to remove this form of U{1) x C/(l) symmetry 
breaking and the HQV lattice structures are still possible 
(indeed the theory is the same as that given for the non- 
chiral spin-triplet superconductors with a Zeeman field, 
but without spin-orbit coupling). Note that if t ^ 
(signaling the existence of the fractional vortex lattice), 
then (|^oP0o02)= (l</'oP02<^o)='-' ^^"^ ^"^y lattic e ge ome- 



try. Consequently, the last two terms of Eq. (6.5 1 ) do 



not play any role in the theory of the fractional vortex 
lattices. 

It is reasonable to ask if there are any other U{l)xU{l) 
symmetry breaking terms that we have neglected in the 
above analysis. Indeed there is one that appears at or- 
der e^: t^l3i'^l{-il)*{(j)l{^ly). This term allows for the 
existence of a fractional vortex lattice phase subject to 
the constraint that t is half a vortex lattice translation 
vectoil^. There are also U{\) x U{\) that appear at or- 
der e^, but these vanish for the same reason as the order 
e term. Consequently, the spin-orbit coupling for the 
hexagonal two-dimensional representations plays essen- 
tially the same role as the Zeeman field. 



VII. OBSERVATION OF THE VORTEX 
LATTICE 



The best way to determine both the vortex lattice 
structure and the vortex type is to observe the magnetic 
field distribution through the small angle neutron scat- 
tering. What we will see in this experiment is the Fourier 
transform /(G) of the screening field of Eq. (5 



hs{v) = ( 



+ ( 



$0 



47r«)C2|0o(r)P 
47r«)C2|0„(r)|2 



(7.1) 



- that is hs{v) = ^Q/(G)exp(zG • r), where G is the 
reciprocal lattice vectors of the vortex lattice in the unit 
of the inverse magnetic length. 

The characteristic feature of the vortex lattice with 
half- vortices in the small angle neutron scattering exper- 
iment is the modulation of the Bragg peaks. The form 
factor of the Bragg peaks would be the /(G) of the last 
paragraph. Using 



/'o(r)P = 5](-l) 



mi+m2-|-mim2g— G /2 



(7.2) 



where G = rniGi -|-m2G2, and Gi's are the basis vector 
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FIG. 4: The contour plots of the screening field in the real (that is, position) space for different types of vortex lattices: (a) 
a lattice of ordinary Abrokosov (full quantum) vortices (b) a HQV lattice (c) a lattice of fractional vortices in the presence of 
Zeeman field. The position is measured in units of magnetic length. Note the halving of the unit cell in going from the HQV 
lattice to the full quantum vortex lattice. When the Zeeman field is added in, the periodicity is that of the full quantum lattice, 
though the vortex lattice unit cell now has additional structure due to appearance of fractional flux at both the corners and 
center of the unit cell. 



of the reciprocal lattice, we obtain the form factor 

/(G) = (-l)"i+™^+'"i"^e^^'/2[(^^ -47rK)|CttP 

+ (^ + 4.^)|Cj,pe'«-]. (7.3) 

This equation implies that the intensity |/(G)p for our 
Bragg peaks does not come out same for all G's. This 
is because for almost all vortex lattice structure (the 
single exception being not very robuts honeycomb lat- 
tice) T is half a vortex translation vector so that we have 
e'*^"^ = — 1 for half of G's and e'*^"^ — 1 for the other 
half. When there is no Zeeman field, e**^ "^ = — 1 peaks 
disappear completely; natural given that magnetic field 
cannot distinguish the spin up and the spin down HQV's 
at all and thus sees the unit lattice vector halved. How- 
ever, when the Zeeman field breaks down the Z2 symme- 
try between the spin up-up pairs and down-down pairs, 
we now see a secondary peak for e**^"^ = — 1 as shown on 
Fig. (§. 

Another promising direction for detecting fractional 
vortex lattice would be to use spin-polarized STM to 
probe the vortex cores. The key point is that the low 
energy quasi-particle spins have opposite polarization in 
the two different HQV's. This is because for half of HQV 
cores, we have A|| = and A|| ^ 0, so that only spin- 
down quasi-particles are gapped. On the other hand, for 
the other half of HQV cores, only spin-up quasi-particles 
are gapped. This spin-polarization of the subgap core 
modes should be readily detected through spin-polarized 
STM. 



VIII. CONCLUSION 

In this paper we explored various possibilities for frac- 
tional vortex lattice structures in spin triplet supercon- 



ductors starting from the most general from of Gibbs 
free energy that is allowed by the symmetry of the or- 
der parameter and that of the lattice symmetries rel- 
evant for three candidate spin triplet superconductors, 
namely single layer ruthenate Sr2Ru04^^ cobaltate 
Na:rCo02 • 2/H20C3 and organicP (TMTSF)2C104. The 
focus of our analysis was on the role of aspects unique 
to triplet superconductors, such as (i) Cooper pair Zee- 
man field, (it) spin-orbit coupling, (iii) screening, and (iv) 
interaction effects in the energetics of the vortex lattice 
structure, (i) The Cooper pair Zeeman field breaks a Z2 
symmetry of the free energy whose presence constrains 
the fractional vortices to contain half integral flux quanta. 
The resulting structure is that of two interlacing lattices 
of vortices containing arbitrary fraction of flux quanta 
that adds up to one flux quanta. Such fractional vortex 
lattices will have interesting field distributions in vor- 
tex lattice unit cell due to internal structures within the 
unit cell, (ii) The effect of spin-orbit coupling is lattice 
symmetry specific. In hexagonal lattices systems such as 
cobaltates Naa:Co02-yH20, spin-orbit coupling has the 
same effect as the Cooper pair Zeeman field, support- 
ing fractional vortex lattices. However, for tetragonal or 
orthorhombic lattices, sufficiently strong spin-orbit cou- 
pling generally favors ordinary Abrikosov vortex lattice 
over HQV's. However, such an effect is relatively mild 
in a dense vortex lattice, provided that the separation 
between the HQV vortices is less that a length set by the 
spin-orbit coupling, (iii) The Meissner screening effec- 
tively generates attraction between two HQV's with op- 
posite winding of the spin phase and weakly destablizes 
the HQV's. (iv) The interaction effects clearly support 
energetic stability of HQV's within the GL theory. The 
interaction effects represented by inhomogeneous (unique 
to triplet superconductors) quartic terms can drive dif- 
ference in effective superfluid stiffness psp < Ps which 
stabilizes HQV's in the London limit. When the above 
effects are put together, all weak coupling theories we 
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examined appears to lie at the point of fine balance be- 
tween ordinary Abrikosov vortex lattice and lattices of 
HQV's. Hence it should be possible to observe transi- 
tions between these structures with small changes of pa- 
rameters. This further motivates experimental search for 
these fractional vortex lattices. We have sketched possi- 
ble routes for such searches using neutron scattering or 
spin polarized STM. 
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APPENDIX A: GINZBURG LANDAU ENERGY: 
FOURTH ORDER TERMS FROM 
WEAK-COUPLING THEORY 

The GL free energy can be determined in the weak- 
coupling limit. In the context of the existence of 1/2 qv 
lattice structures, the result for the fourth order terms 
in the free energy turns out to be highly relevant. As 
shown here, this reveals that weak-coupling theories sit 
at a point in which the 1/2 qv and the full qv lattices 
are degenerate. This indicates that interactions beyond 
the weak-coupling limit are essential to determining the 
which lattice structure actually appears (screening plays 
a role here as well as shown earlier). 

The portion of the free energy we calculate here is given 
in Eq. 3.6 



(Al) 



The weak-coupling limit (without spin-orbit coupling) 
yields (this follows from Ref. 



/(l(x(|d(k)n + (q^(k)) 



(A2) 



where q(k) = id(k) x d*(k), (/i(k)) means average /i(k) 
over all k on the Fermi surface, and the proportional- 
ity constant can be found but it is not important for 
our considerations. When q is non-zero, then the super- 
conducting state is called non-unitary. In weak-coupling 
theories, non-unitary states cost energy and typically do 
not appear. Using the gap structure of Eq. 2.2, we find 



/•(4) 
•> hom 



(x(|/(k)|4)[(|Anr 
= 2(|/(k)r)(|Anr 



\Au\r+ 

\Aii\'). 



\A 



TTl 



|AxJ^)] 
(A3) 



This implies that f32 = — 2/3i, independent of the shape 
of the Fermi surface. The lack of interaction between the 



two components of the gap function leads to the degener- 
acy between the 1/2 qv and the full qv lattice structures. 



APPENDIX B: RUTHENATE - THE LANDAU 
LEVEL MIXING 

We show here how we can have the Landau level mix- 
ing in a chiral triplet superconductor. The case we are 
considering here is in the weak pairing regime and has 
tetragonal crystalline symmetry and a cylindrical Fermi 
surface. Let us consider again the linearized GL equa- 
tion: 



A.,_ 



K / 1 



2n+n_ 
ni 



1 + 2n+n_ 



A.+ 
A.,_ 



(Bl) 

where s ^TTjU- (Note that, though otherwise same as 



Eq. (5.201, we now ignore the energy splitting between 



the ± chiralities and set Ki = K2 = K.) This matrix 
equation as a solution in the form 



(B2) 



where 5 — a/3 — ^/2. 

When we ignore the Zeeman field, much of the vortex 
lattice energetics of the lowest Landau level case remains 
valid with the Landau level mixing. For instance, the 
two main formulas of Section V. A, Eqs.(5.12l, 



(/) = ^a{Cl + Cl) + /3(C4 + CD + P^CM, (B3) 



and (5.141) 



(/) = - 



8^ 2/3 + /33 



(B4) 



remains valid, mainly due to |A| c>c (-ffc2 — HY/"^. This 
means we can still calculate h^, solely from quadratic 
terms. For quadratic terms, we simply have two copies 
(for s —X\ and ||) of what was obtained for the case of 
d — (fcj, + iky)2 by one of u^, we can use the formula 
for hs for that case: 



= ^^[Cf^{(l-3,5/V2 + 2(5^)1001' 

J/\/2)|0i|2-|-j2|^2|2} 



$0 

+ (2J2 

+ Cf^{(l- 3(5/72 + 2(52) 1 



+ (2(5^- VV2)|(^i|2 + 52|^^|2|]^ 

The spatial average of this equation is still proportional 
to (C^^ + C 
still stands: 



to {C^^ + Cf^) just like Eq. (|5^. Also, a oc (iJc2 - H) 



(2(5^ 



$0 

'6/V2){\4>i\' 



5'm% (B6) 



11 



However, the formula for 2(3 + are much more com- 
plicated here, especially when we include all terms of 



Eqs. (3.14) and (3.15) for these coefficients. For sake 



of convenience, instead of directly writing down /3 and 
/33, we will list h"^ (terms that are proportional to K^), 



fhora (terms involving coefficients of Eq.(3A4)), and _ 
(terms involving coefficients of Eq. ( 3.15[ )); to obtain 
2/3 + /33 , we can use the relation 



(4) 



Stt' 



(B7) 



The following is the full listing of h'^/Sn, fjj^J„^, and /, 
(note that we have set (3^ of Eq. (3.14) to be zero): 



(4) 



[{l-3S/V2 + 26Y{\M') 



J horn 



2(2,52 - V\/2)(l - i5l^/2 + 25-'){\U^\cj)i\^) 
25\\-?,5/^/2 + 25^){\cp^\'\cp2?) 
{25^ ~5/V2fm^) 

25^25^ - 5/^)m^\^2?) + 
{\^?>5/V2 + 26^fm''\~M^) 

2(252 _ ^/72)(1 - 35/V2 + 2,52)(|0^|2|^^|2^ 

25\l-i5/V2 + 25^){\U^\~^2?) 

[25^ -5lV2fmX^^\^) 

26\25^-5/^/2){\<j,,\X^2?) 

5\m^\~^2?)]. (B8) 

m\M')+s\\ct^2\'))+2p[5wu^)m') 



25^l3'^{\U''\h?) 



(B9) 



and 



p(4) 



2i 



-(352(|0o|2|^2|2) 



35'(|0On03p) 



2^4(101^01 

3<5M02n02 



4/IJ, |2|I |2\ 



4/,. ,2,. ,2^ „3^4^|^^|2|^^| ^ 

354(102^103^). (BIO) 



APPENDIX C: CORRELATION FUNCTIONS 



In calculating (/f^^), note 



(£>/) . {Dgf = ^[(n+/)(n+5)* + (n_/)(n_g)*] (ci) 



and 



n+0„ = \/Tr+T0„+i 

Together with partial integration 



(C2) 



((n+0„)0*j0p0*) = (0„(n_0,„)*0p0*) 

- (0nCn(n+0p)0q) 

+ (0„0:;,0p(n_0,)*), (C3) 



these equation gives 



^(0S0o(r>0o)-(D0o)*> = (|<^oPl<^oP>-(|0on<^i|') 



1 



6S02(r>0o) • (D02)*> - ~ {\M'm') 



^02(r>02) • (r>02)*) = 2(|0o|l0l|^) + (100110212) 
+ 2(|0in0i|2) + (|0i|2|02|2) 
- 3(|0on03|2)-3(|0in03|2) 

+ 3(|02n02|2)-3(|02n03p). 

(C4) 

These can be evaluated using 



(l'^0|2)2 



Y^Ll{klj2)Ll{klj2)e-'^'J\ 

r,s 

= Ll{klj2)Ll{klj2)e''''r^'^ cosikrs ■ t) 

(C5) 
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